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The interaction of exciton and charge transfer (CT) states plays a central role in photo-induced CT 
processes in chemistry, biology and physics. In this work, we use a combination of two-dimensional 
electronic spectroscopy (2D-ES), pump-probe measurements and quantum chemistry to investigate 
the ultrafast CT dynamics in a lutetium bisphthalocyanine dimer in different oxidation states. It 
is found that in the anionic form, the combination of strong CT-exciton interaction and electronic 
asymmetry induced by a counter-ion enables CT between the two macrocycles of the complex on a 
30 fs timescale. Following optical excitation, a chain of electron and hole transfer steps gives rise 
to characteristic cross-peak dynamics in the electronic 2D spectra, and we monitor how the excited 
state charge density ultimately localizes on the macrocycle closest to the counter-ion within 100 fs. 
A comparison with the dynamics in the radical species further elucidates how CT states modulate 
the electronic structure and tune fs-reaction dynamics. Our experiments demonstrate the unique 
capability of 2D-ES in combination with other methods to decipher ultrafast CT dynamics. 



I. INTRODUCTION 

Charge transfer (CT) is a ubiquitous process in na- 
ture representing the initial step in many chemical re- 
actions and bioenergetic pathways.^ Photo-induced CT 
is intimately related to harnessing and conversion of ra- 
diative energy in photovoltaics^^*, photosynthesis^"^ and 
catalysis*. CT states, representing an excitation in a 
molecular aggregate where the electron and hole are lo- 
cated on different chromophores, are characterized by 
low transition dipole moments and strong sensitivity to 
environmental perturbations.^'^" Due to their low tran- 
sition dipole moments, CT states participate in photo- 
induced reactions via the interaction with optically al- 
lowed (bright) states. If the interaction with the CT 
states is weak, the resonance interaction between locally 
excited (LE) states leads to excited states which are de- 
localized over the chromophores in the aggregate (exci- 
ton states). For a more pronounced CT interaction, the 
mixing of excitonic and CT states leads to (delocalized) 
eigenstates exhibiting various degrees of charge separa- 
tion. After photo-excitation, interaction with the nuclear 
degrees of freedom enables relaxation within the manifold 
of excited states, ultimately leading to population of the 
CT states. 

The influence of the CT states on the electronic struc- 
ture depends on the relative magnitude of resonance- 
(Coulombic interaction between LE states) and CT- 
coupling (between LE and CT states, depending on the 
wave function overlap^^). This elementary interplay can 
be studied in model systems where the different cou- 



plings can be controlled via the inter-pigment distance, 
and where the energies of the CT states can be altered 
electrochemically. One such class of molecules is bisph- 
thalocyanines, which comprise two macrocyclic phthalo- 
cyanine moieties held together at a flxed distance by a 
rare earth cation. ^^'^^ The close proximity of the two 
rings (2.8 A for lutetium) gives rise to delocalized charge 
distributions and results in strong CT coupling due to 
significant wavefunction overlap. ^^ The resulting high 
degree of charge carrier mobility readily explains their 
semiconductivity^^"'^^ and electrochromic behavior^^'^*. 

In its electronic ground state, lutetium bisphthalocya- 
nine ([LuPc2]) is a stable radical ([LuPc2]*) with an inter- 
mediate oxidation state of —1.5 on each ring^^. [LuPc2]* 
cannot stabilize inter-ring CT states, and the electronic 
structure can be understood via resonance interaction be- 
tween the LE states. The electronic structure will be sim- 
ilar to that of an H-type dimer, which gives rise to a single 
transition in the linear absorption spectrum (Fig. 1). On 
the other hand, the anion ([LuPc2]~) exhibits CT states 
in the same energy range as the LE states. ^^ The ad- 
ditional CT coupling drastically changes the electronic 
structure and gives rise to two bands in the linear ab- 
sorption spectrum (Fig. 1). Therefore, the comparison 
of [LuPc2] in different oxidation states provides a perfect 
opportunity to pinpoint the infiuence of the CT states on 
the electronic structure and dynamics. 

The structural and electronic similarities favor 
[LuPc2] ^ as a biomimetic model system for studies of the 
special pair in the photosynthetic reaction center (RC)"''. 
However, for isolated [LuPc2]~, the CT states represent- 



ing forward and backward CT between the rings are en- 
ergetically degenerate, and the electronic structure can 
be described by charge resonance states. ^^ The presence 
of an environmental asymmetry breaks the energetic de- 
generacy of the CT states, lowering the energy of the 
CT state on one side of the complex and enables a net 
transfer of charge after photo-excitation. In the RC, such 
asymmetry arises from differences in the electronic cou- 
pling among pigments due to slightly different spacings 
between the cofactor units^°, distinct dielectric environ- 
ments due to intra-protein electric fields^^"^"^, and differ- 
ent hydrogen bonding pattern around the central units^^. 
In other redox systems, the asymmetry may be provided 
by counter-ions. The formation of ion associates can have 
profound influence on the electronic structure and the 
path of chemical reactions. ^'"'"■^^ For weak CT coupling 
(non-adiabatic CT), the role of the counter-ion has been 
studied, and it has been shown that the counter-ion can 
actively or passively control the rate of CT depending on 
the strength of association. ^^"^^ 

Conventionally, CT has been studied in systems where 
the CT coupling is larger than the resonance coupling, 
implying that CT proceeds without interference from en- 
ergy transfer processes. However, in systems like the 
natural RCs or in novel architecture for artificial light 
harvesting, the two couplings can be of equal magni- 
tude and both energy transfer and CT will take place 
simultaneously. The interplay of resonance- and CT- 
coupling in these systems tunes the excited state evo- 
lution and the dynamics of the charge separation and 
energy transfer processes. Such excited state dynamics 
gives rise to complicated spectral signatures, making a 
detailed investigation of the underlying mechanisms with 
conventional techniques difficult. In this work we em- 
ploy two-dimensional- (2D) and ultra-broadband pump- 
probe-spectroscopy to disentangle the excited state dy- 
namics in a metal bridged dimer ([LuPc2]), where the 
resonance- and CT-couplings are of equal magnitude. 
The infiuence of CT states and CT coupling on the 
excited state structure and dynamics is further stud- 
ied via electrochemical tuning of the CT state energies. 
In addition, we investigate how specific details of ion- 
pairing enable ultrafast charge transfer across the com- 
plex. The combination of various approaches, both the- 
oretical (quantum chemistry and density matrix propa- 
gation) and experimental (absorption, pump-probe and 
2D spectroscopy) , gives a detailed picture of the ultrafast 
CT dynamics. 



II. EXPERIMENTAL METHODS 

Lutetium(III)bisphthalocyanine-tetrabutyl- 
ammonium salt ([LuPc2]~TBA+), and its radical 
form [LuPc2]* were synthesized via a base catalyzed 
cyclotetramerization reaction of phthalonitrile with 
lutetiumacetate in the bulk-phase at elevated tem- 
perature. The synthesis was carried out using known 
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FIG. 1: Linear absorption of [LuPc2]^TBA+ in benzonitrile 
(blue line) and [LuPc2]* in toluene (red line) shown together 
with the pulse spectrum used in the 2D experiments (grey). 
The inset shows the structure of [LuPc2]^ and the position of 
the TBA+ counter-ion. 



procedures from the literature. ^■^'■^° Both the synthetic 
preparation and the characterization of the complex are 
described in detail in the Supporting Information.^^ 

The setup and methodology used for the 2D experi- 
ments has been described in detail previously.^^ In the 
experiments reported here, a home built NOPA oper- 
ated at 200 kHz repetition rate generated pulses centered 
at 15000 cm-i with a bandwidth of 2500 cm"! FWHM 
(Fig. 1). The FWHM of the intensity autocorrelation 
was 13 fs (9.2 fs pulse duration), while the correspond- 
ing width assuming a flat phase over the pulse spectrum 
was 10 fs. Two phase stable pulse pairs were generated 
by diffracting the NOPA pulses through a transmission 
grating. In the 2D experiment, the delay between the 
flrst two pulses, ti, was scanned with interferometric ac- 
curacy, while the delay between the second and third 
interaction, ^2, was held constant. The signal radiated in 
the phasematching direction was overlapped with the lo- 
cal oscillator (LO) , and the signal field was recovered by 
spectral interferometry. The 2D spectrum as a function 
of vi and V3 was obtained by Fourier transforming over 
ti and addition of rephasing and non-rephasing signal 
contributions. The absolute phase of each 2D spectrum 
was determined by adjusting the LO delay and phase to 
optimize the overlap of the projection of the real part 
of the 2D spectrum with the spectrally resolved pump- 
probe signal. ^^ [LuPc2]~TBA"'" dissolved in benzonitrile 
was used in a wire-guided fiow jet giving a film thickness 
of 230 /im. The structure and linear absorption spectrum 
of [LuPc2]~TBA+ are shown in Fig. 1 together with the 
NOPA spectrum used in the experiments. 

Broadband pump-probe measurements were carried 
out in a setup described in detail previously. ^'^ To ex- 
tend the probe range to the NIR region, a home built 
OPA operating at 8300 cm~^ was used to generate white- 
light from 9000-25000 cm-^.^'^'^ To cover the bands in 



the UV region (25000-32000 cm-i), white-light was gen- 
erated by focusing the output from the Ti:Sapphire laser 
into a CaF2 plate. ■^'^ Pump pulses of approximately 40 fs 
duration were used to selectively excite the two bands at 
14200 and 16200 cm~^ giving a time resolution of 50 fs. 



III. RESULTS 
A. Experimental Results 

The interplay of exciton and CT states and their influ- 
ence on the excited state dynamics in [LuPc2]^ have been 
studied previously with pump-probe- and two-color pho- 
ton echo-spectroscopy.'^^''^^ The pump-probe measure- 
ments revealed that the transient spectra were indepen- 
dent of which of the two bands was pumped, indicat- 
ing rapid relaxation in the excited state. To resolve the 
excited state dynamics we turned to 2D spectroscopy."^^ 
By correlating the optical coherence evolution between 
the first and third time interval in a four-wave mixing 
sequence, 2D spectroscopy maps the system dynamics 
onto two dimensions, thereby minimizing spectral con- 
gestion and provides a very detailed picture of the sys- 
tem's evolution.'^* 

Fig. 2 shows the real part of the 2D spectrum of 
[LuPc2]^TBA+ in benzonitrile for t2 = 0, 15, 45, and 
400 fs. The waiting time t2 denotes the delay between 
the first and the second pulse pair, and corresponds to 
the probe delay in a pump-probe experiment. All spec- 
tra have been normalized to their respective maxima, and 
contour lines are drawn in 5 % increments. In Fig. 2, pos- 
itive signals correspond to stimulated emission (SE) and 
ground state bleach (GSB), while excited state absorp- 
tion (ESA) results in negative signals. The 2D spectrum 
at ^2 = 15 fs (outside pulse overlap) shows two diagonal 
peaks corresponding to the two bands in the linear ab- 
sorption spectrum. A closer inspection of the low energy 
diagonal peak reveals that it is more accurately described 
as two separate peaks at v^ ~ 14200 and 14800 cm^^, 
respectively. Above the high energy diagonal peak we 
observe small amounts of ESA (blue dashed lines). In 
addition to the diagonal peaks and ESA contributions, 
one cross-peak at vi, v^ = 14200, 16200 cm~^ is clearly 
visible. This cross-peak shows that the transitions re- 
sponsible for the main bands in the linear absorption 
spectrum share a common ground state. Surprisingly, 
there is no corresponding cross-peak below the diagonal. 
For longer i2 we find that the cross-peak below the diag- 
onal recovers. However, the recovery of the cross-peak is 
not uniform. The 2D spectrum at i2 = 45 fs shows that 
the cross-peak at v^ = 14800 cm^^ grows much faster 
than the one at 14200 cm~^. For longer ^2, the cross- 
peak below the diagonal becomes more symmetric, but 
its amplitude remains weaker than the amplitude of the 
opposite cross-peak for all t2- 

The 2D spectra in Fig. 2 do not provide an accurate 
picture of the dynamics in the system because of the 



normalization. The variation of the amplitudes of the 
different peaks as a function of ^2 provides information 
on population dynamics of the system which cannot be 
obtained from the normalized spectra. To analyze the 
kinetics, we show in Fig. 3 the time evolution of the vol- 
ume of a 30 X 30 cm^^ box for the diagonal peaks at 14200 
(dpi), 14800 (dp2) and 16200 cm"! (dp3), and the cross- 
peaks wi, W3 = 16200, 14300 cm^i (cpl), wi, W3 = 16200, 
14800 cm-i (cp2), and Vi, wg = 14200, 16200 cm^^ (cp3). 
The positions of these points are also shown in Fig. 2 for 
^2 = 45 fs. The kinetics of the diagonal peaks all dis- 
play a fast decay to about 50 % of the value at <2 = 0. 
The amplitude of dp3 decays on a 30 fs timescale, while 
dpi and dp2 show an almost constant amplitude outside 
pulse overlap (^2 > 15 fs). The cross-peaks exhibit more 
complicated behavior. While cp3 behaves qualitatively in 
the same way as dp3, the cross-peaks below the diagonal 
show distinct dynamics. The amplitude of cp2 rises on a 
20 fs timescale followed by a quick decay. After about 50 
fs, the rise is complete and the signal shows a slow oscil- 
lation around a stationary value. The amplitude of cpl 
shows a delayed rise, reaching its peak value after about 
100 fs. After reaching their final values, cpl and cp2 
both oscillate with a frequency of 160 cm~^. The same 
modulation is found in pump-probe"^^, two-color photon 
echo^^. We can assign this mode, based on our quan- 
tum chemical calculations, to a modulation of the Lu-N 
distance. ^^ 

The appearance of the 2D spectrum for ^2 > 30 fs 
deviates from what we expect for a coupled dimer, in- 
sofar as it does not display a characteristic four peak 
pattern"^^. The absence of a cross-peak below the diago- 
nal gives rise to a pronounced asymmetry of the spectra. 
This effect clearly remains outside pulse overlap, but di- 
minishes as the population relaxes from the upper band 
(decay of dp3). This points to ESA from the upper band 
as the origin of the cancellation of the cross-peak below 
the diagonal. As the ESA signal decays, the underly- 
ing GSB contribution becomes visible. Such relaxation 
should give rise to a SE signal in the region of the cross- 
peak below the diagonal (relaxation cross-peak) . For this 
reason, we expect this cross-peak to be stronger than the 
corresponding one above the diagonal (cp3) when relax- 
ation from the upper band has been completed. However, 
we find that the cross-peak below the diagonal remains 
weaker for all ^2- This observation could in principle be 
explained by loss of the SE signal due to relaxation to a 
state outside the spectral window or into a state with low 
transition dipole moment. However, we see no indication 
of such relaxation in the kinetics of either dpl+dp2 or 
cpl+cp2. If the population remains in the lower band, 
the positive SE signal could be hidden by a negative ESA 
signal overlapping with the positive GSB of the low en- 
ergy band. While this would explain the 2D kinetics, it 
cannot explain the reported lack of fluorescence from the 
low energy band"^^. 

It thus becomes apparent that the 2D experiment, de- 
spite the 4500 cm~^ probe range, is not able to follow all 
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FIG. 2: 2D spectra of [LuPc2]^TBA+ at different waiting times (t2) indicated in the Figure. The top row shows the experimental 
results, while the middle and bottom rows show the simulated total spectra and ground state bleach (GSB) plus stimulated 
emission (SE) contributions, respectively. All spectra have been normalized to the maximum of the total signal, and contour 
lines are drawn in 5 % increments starting at 5 %. Positive signal contributions (GSB+SE) are shown in red while negative 
(ESA) are shown in blue dash. The points used for the kinetics in Fig. 3 are shown for t2 = 45 fs. 



relaxation processes in the system and cannot provide on 
its own the basis for a complete relaxation model. To be 
able to probe the response of (possible) low energy states 
in the NIR which are not covered in our 2D experiment, 
we turned to broadband pump-probe spectroscopy. Fig. 4 
shows the pump-probe spectra in the NIR region for a few 
selected delays. At ^2 — 100 fs, a broad ESA signal covers 
the entire region from 9000-13000 cm~^. For longer t2, 
the ESA around 12500 cm~^ (blue arrow) seems to de- 
cay significantly faster than the ESA below 11000 cm~^ 
(red arrow). The single wavelength fits shown in the in- 
set reveal that the ESA at 10870 cm~^ decays with a 
3.8 ps exponential component, which is equal to the life- 
time of the lowest excited state obtained from a global 
fit.^^ On the other hand, the signal at 12500 cm^^ shows 
an additional 70-100 fs decay component. The kinetics 
in this spectral region reflects the sum of ESA and SE 
contributions. The additional fast decay of the signal 
at 12500 cm^^ can thus be interpreted as the delayed 



rise of SE from an almost dark state 
in the NIR region as discussed above. 
12500 cm~^ is about 20 times weaker 
20 mOD) than the main GSB peak at 
low transition dipole moment of this 
explains the lack of fluorescence in this 
ps lifetime, we estimate a fluorescence 
the 10~^ range, beyond the detection 
and previous experiments^^. 



with a transition 
The SE signal at 
('-^l compared to 
14200 cm-i. The 
transition readily 
system. For a 3.8 
quantum yield in 
sensitivity of our 



The discussion of 2D and pump-probe spectra above 
provides a rough overview of electronic structure and ul- 
trafast dynamics in [LuPc2]^TBA+. The states initially 
prepared by pumping the upper band decay to the low 
energy bright band on a 30 fs timescale. The presence of 
two diagonal peaks, and the distinct relaxation dynam- 
ics associated with cpl and cp2 points to two separate 
electronic states in this band (Fig. 3). The population of 
these two states flows on a 100 fs timescale into an almost 
dark state with a transition in the NIR range (Fig. 4). 
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FIG. 3: Kinetics of the different points in the 2D spectra indi- 
cated in Fig. 2. (a) Diagonal peak kinetics for dpi (black tri- 
angles), dp2 (red circles) and dp3 (blue squares, upper panel) 
together with the simulations (solid lines in the same color). 
The different curves have their own y-axis as indicated by 
the colors, (b) Cross-peak kinetics for cpl (black triangles), 
cp2 (red circles) and cp3 (blue squares, upper panel) together 
with the simulated curves (solid lines). 
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FIG. 4: Pump-probe spectra of [LuPc2]^TBA^ in benzoni- 
trile with a focus on the NIR region for f2 ~ 0.1, 0.6, 1.5, 3, 
and 10 ps (from red to blue). The inset shows the (normal- 
ized) kinetics at the frequencies indicated by the arrows in 
the main Figure, i.e. at 10870 cm~^ (blue circles) and 12500 
cm~^ (red circles). The solid lines are the corresponding fits. 
In pump-probe spectroscopy, the standard sign convention 
(i.e. AOD) is opposite to that in 2D, e.g. ESA gives rise to 
positive signal contributions. 
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FIG. 5: Pump-probe spectra of [LuPc2]* in toluene for f2 = 
0.05, 0.2, 0.5, 1, 3, and 5 ps (from bottom to top). The spectra 
have been shifted for ease of viewing. Solid lines are fitted re- 
sults based on a model discussed in the text. Negative signals 
correspond to SE/GSB in accordance with the convention in 
pump-probe. 



To elucidate the influence of the CT states on the elec- 
tronic structure in [LuPc2]; we carried out pump-probe 
experiments on [LuPc2]*. Fig. 5 shows pump-probe spec- 
tra at selected delays together with the fit based on a 
kinetic model including the upper (optically allowed) ex- 
citon state, the lower (optically forbidden) exciton state, 
and a hot ground state. From the fit we conclude that 
the upper exciton level decays to the lower one with a 400 
fs time-constant. From this level, there is an ESA transi- 
tion to a doubly excited state. This doubly excited state 
is shifted by —550 cm~^ as compared to twice the LE en- 
ergy. The red-shifted ESA from the dark state leads to 
a seemingly faster decay of the signal on the red side of 
the spectrum as the population fiows into the low energy 
exciton state. Finally, the lower exciton level decays into 
a hot ground state with a 2.5 ps time-constant, where 
subsequent cooling takes place on a timescale of 15 ps. 



B. Electronic Structure of [LuPc2] 

To be able to interpret the experimental results on 
[LuPc2] and specifically to elucidate the role of the CT 
states, it is necessary to develop a model of the elec- 
tronic structure enabling simulations of the linear- and 
nonlinear-spectra. Such a representation is provided 
by the Frenkel exciton Hamiltonian with LE and CT 
states. ^*^ This model provides a descriptive picture of the 
interaction of optical excitations with CT states, which is 
useful for generalization and interpretation of the results. 
Furthermore, it allows us to include doubly excited states 
needed to account for ESA, which is extremely important 
for the interpretation of the 2D experiments. 

As a starting point we include the Q^.- and Qy- 
excitations on each monomer (A and B). [LuPc2] exhibits 
Did symmetry^^, e. g., the two rings are twisted by 45 
degrees with respect to each other (Fig. 1). The tran- 
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FIG. 6: Diagrams illustrating selected singly excited states 
contributing to the model. The full set of states is obtained 
by permutations of A — > B. Solid circles represent electrons, 
hollow circles holes. Black and red color refers to excitation 
of X- and y-orbitals, respectively. The couplings between the 
states are indicated by Jex and Jet- 



sitions on the different monomers interact via resonance 
coupling, and due to the overlap of the wave functions, 
the excitations on the different rings are also coupled via 
electron transfer into the LUMO orbitals. This allows 
us to construct 4 distinguishable CT states. In our ter- 
minology, C^ is the state where the excited electron is 
transferred from the excited x- or y-orbital on monomer 
B to the excited x-orbital on monomer A. A subset of 
the states used to construct the basis states is shown in 
Fig. 6. Using the configuration diagrams we can write 
the Hamiltonian for the one-exciton manifold as 



Hs = 
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where J^^ is the resonance- and Jd 
Taking the rings as equivalent, Qx 



(1) 
is the CT-coupling. 

Qy, and assuming 



that C^ > Q^ , the present model recovers the previously 
proposed electronic structure of [LuPc2]~^^''*^''*^. In this 
model, the band at 14200 cm~^ has predominant charge 
resonance character, while the band at 16200 cm~^ has 
a dominant exciton character. In addition, this model 
predicts one dark transition in the NIR region.^* This 
exciton model agrees very well with our TD/DFT cal- 
culations on [LuPc2]^ discussed in Appendix A. These 
calculations yielded three doubly degenerate electronic 
states with vertical transition energies of 11682, 15848 
and 18657 cm"'^, and oscillator strengths of 0.00, 0.45 
and 1.03, respectively. 

Both the exciton model in the form discussed above 
and the direct TD/DFT calculations predict two bright 
states in the relevant spectral range, i. e. the two main 



FIG. 7: Optimized ground state geometries 
BHLYP/SV(P) level of theory. (a) [LuPca]" 
[LuPc2]"TBA+. 
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absorption bands. However, none of the models are able 
to explain the double peak structure of the low energy 
band or the distinct dynamics of the cross-peaks below 
the diagonal in the 2D spectra. To find an explanation 
for these additional features found in the experiment, we 
investigated the role of the counter-ion (TBA"*") present 
in our experiments. NMR studies on a series of bisph- 
thalocyanines have shown a considerable amount of ion 
association in solution. '^^'^'^ These experiments concluded 
that the counter-ion is situated on top of one of the Pc 
macrocycles as indicated in Fig. 1. This is similar to the 
arrangement in the crystal structure^^, which we took 
as a starting point for the calculations of the optimized 
ground state geometries and vertical transition energies 
of the [LuPc2]^TBA+ complex. The calculations show 
that the interaction with the counter-ion leads to a split- 
ting of the doubly degenerate states found for [LuPc2]~, 
giving rise to 6 transitions in the relevant spectral range. 
The energies, oscillator strengths and contributing or- 
bitals for [LuPc2]~TBA+ as derived from quantum chem- 
istry are shown in Table I. 

The association with the counter-ion leads to signifi- 
cant changes in the electronic structure of the complex 
and gives rise to additional allowed electronic transitions. 
To understand how these changes can be interpreted in 
terms of the energies and couplings between the states 
in the exciton model. Fig. 7 shows the optimized ground 
state geometry for [LuPc2]^ and [LuPc2]^TBA+. For 
[LuPc2]~, the di-anionic rings (Pc^~) display a symmet- 
ric off-planarity distortion due to electrostatic interaction 
with the central metal ion. In the presence of the counter- 
ion, the off-planarity distortion is decreased significantly 
along the axis in contact with the alkyl chains of the 
TBA+ molecule. From an exciton model point of view, 
this implies that the excitation along this axis, here Qj;, 
on the ring closest to the counter-ion should be lowered 
in energy. Similarly, the state corresponding to CT into 
the same orbital should decrease in energy as well. 

The simulations of the linear absorption spectra of 
[LuPc2]* and [LuPc2]~TBA"'" based on Eq. 1 are shown 
in Fig. 8. For [LuPc2]^TBA+ we also show the stick- 
spectrum obtained after direct diagonalization of Hg . De- 
noting the ring closest to the (positive) counter-ion A, 
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TABLE I: Transition energies and oscillator strengths for the first 6 electronic states in [LuPc2] TEA"'" at the BHLYP/SV(P) 
level of theory. 



we find that lowering the energy of C^ (electron trans- 
fer towards the counter-ion) is able to mimic the general 
features of the electronic structure found in the quantum 
chemical calculations. For Fig. 8 we used Q^ — Qy — 



Q^ = Q 



13500 cm-^ C^ 



C' 



C^ = 15150 



cm~^, C^ = 11600 cm^^, Jex — 1395 cm^^, and J^ ~ 
725.5 cm^^. Based on the quantum chemical calcula- 
tions, we expect that the interaction with the counter- 
ion should lead to a decrease of the energy of Q^ as well. 
However, changing the energy of Q^ by —500 cm^^ only 
leads to minor changes in the linear absorption spectrum 
and we therefore kept all Q-transitions degenerate for 
the remainder of the paper. To avoid confusion with the 
states obtained from quantum chemistry, we will denote 
the states of the exciton model with e„. States ei,e2, and 
63 are located in the NIR and have very low or zero tran- 
sition dipole moments. States 64 and 65 are located in the 
low energy band seen in linear absorption and are sepa- 
rated by ~500 cm"'^ (Fig. 8). This double peak structure 
in the low energy bright band is essential to reproduce 
the observations in the 2D spectrum. State eg is dark and 
is located in between the two bands, while states 67 and 
eg are the origin of the upper band in linear absorption. 

With the large number of states and couplings, a 
fit to the linear absorption spectrum is by no means 
unique. An important test of the validity of the model 
and the parameters is to apply the model derived for the 
[LuPc2]-TBA+ to [LuPca]*. The radical itself is non- 
polar and can therefore not stabilize the CT states. We 
modeled the spectrum of [LuPc2]* by raising the ener- 
gies of the CT states to a spectral position where their 
influence on the Q transitions is negligible. In addition, 
we used that the Qx/Q^-transitions are degenerate in the 
absence of the counter-ion. The successful prediction of 
the transition frequency, width of the main peak, and vi- 
brational sideband for [LuPc2]* evident from Fig. 8 shows 
that we have correctly estimated the spectral deirsities, 
resonance coupling and Q-transition energies in [LuPc2]. 

The quantum chemical calculations and the exciton 
model represent two different descriptions of the elec- 
tronic structure of [LuPc2] . The quantum chemical calcu- 
lations adopt a supra-molecular perspective and include 
resonance- and CT-couplings between the two rings im- 
plicitly. The exciton model is based on distinguishable 
molecular excitations with free parameters (energies and 




V (1000 cm-1) 

FIG. 8: Linear absorption spectra of [LuPc2]~TBA+(blue cir- 
cles) and [LuPc2]* (red squares). The solid curves show the 
simulations based on the model discussed in the text. For 
[LuPc2]~TBA^, the stick-spectrum using a single realization 
of disorder is also shown. The four peaks correspond to the 
transitions to the states 64 (14273 cm~^) 65 (14778 cm~^), 67 
(16153 cm~^), and eg (16349 cm^^) in the exciton model. 



couplings) and a minimum number of interactions. Using 
the language of the exciton model, the quantum chemical 
calculations iirclude more states and couplings and will 
provide a more realistic molecular view. For instance, the 
quantum chemical calculations account for hole transfer 
and polarization of [LuPc2]~TBA"'" in the ground state 
not included in the exciton model. However, the quan- 
tum chemical calculations do not provide any insight into 
the dynamics of the system or provide any information 
about doubly excited states, and therefore the combina- 
tion of both models is needed for the interpretation of the 
experimental results. The states of the two models can- 
not be linked one-to-one, but a comparison of the excited 
state manifolds'^^ suggests that one can well connect the 
dynamics of states 64, 65 to (73, 94 and 67,68 to q^ and qq. 



C. Excited State Dynamics in [LuPc2]"TBA+ 

The preceding paragraph outlined a realistic descrip- 
tioir of the electronic structure of [LuPc2]^TBA+. To 



interpret the observed excited state dynamics we car- 
ried out simulations of the 2D spectra based on the CT- 
Hamiltonian (Eq. 1) and time-dependent Redfield theory 
in the Markov approximation (see Appendix B). The sim- 
ulated 2D spectra and kinetics are shown in the middle 
and bottom row in Fig. 2 and Fig. 3, respectively. 

The simulations are in very good agreement with the 
experimental results, and reproduce the shape and evo- 
lution of the 2D spectra. At ^2 = 0, there is a clear 
asymmetry, where the cross-peak below the diagonal is 
significantly weaker than the one above the diagonal. It 
follows from the comparison of the total simulated signal 
with the GSB+SE contribution (Fig. 2) that the miss- 
ing cross-peak is due to overlap of ESA from popula- 
tion in the upper band. The simulations also reproduce 
the double peak structure of the low energy band, and 
we can assign this feature to the eigenstates 64 and 65 
of the Hamiltonian Hg. The splitting between dpi and 
dp2 is close to the frequency of one of the vibrational 
modes found in the quantum chemical calculations (725 
cm^^)^^. If dp2 was a vibrational sideband to dpi, this 
mode would need to have a Huang-Rhys factor on the 
order of 1 to match the observed amplitude. This would 
give rise to a progression of peaks not seen in the experi- 
ments. Furthermore, within the Condon approximation, 
all vibrational transitions involving the same mode would 
have the same direction of the transition dipole moment 
and we would thus expect to see clear cross-peaks be- 
tween them.**^ The lack of a cross-peak between dpi and 
dp2 is readily reproduced in our model (Fig. 2). This can 
can be understood from the perpendicular orientations of 
the transition dipole moments of e4 and e5, which sup- 
presses cross-peaks by a factor of three in an experiment 
with all parallel polarizations.'*^ 

After excitation, coupling to the bath drives relaxation 
between energy levels in the one-exciton manifold lead- 
ing to dynamics in the 2D spectrum. Fig. 9 shows the 
simulated evolution of the populations after excitation 
of state 67 and eg corresponding to the band at 16200 
cm^*. The initially prepared states decay and populate 
state 65 within 30 fs while the population of state 64 rises 
somewhat slower. This is particularly obvious for initial 
excitation of state ey (Fig. 9(b)), where a clear separation 
of the rise of the population of 64 and 65 can be observed. 
By inspecting the relaxation rates, we can conclude that 
relaxation from eg — >■ 65 and 67 — >■ 65 (Redfield rates'^* 
of 18 and 11 fs) is about 3 times faster than relaxation 
to state 64 (64 and 75 fs). This indicates that the faster 
rise of cp2 can be traced to faster population relaxation 
into state 65. To investigate if this relaxation path also 
leads to clear signatures in the 2D spectra, we show in 
Fig. 10 the simulated kinetics with and without the ESA 
contribution. Inspecting the kinetics, we find that the 
cross-peaks indeed show a fast rise reflecting SE after 
population transfer from the upper band, and cp2 rises 
somewhat faster than cpl. However, the SE contribution 
in cp2 does not acquire significant amplitude due to fast 
relaxation from state 65 — > 64. Although the ESA contri- 



bution strongly reshapes the kinetics in the 2D spectra, 
we can conclude that the dominating relaxation path- 
way from the upper band involves relaxation via state 65 
further on to 64. 

Turning to the diagonal peaks, we find that the SE 
signal decays on a 250 fs timescale. Interestingly, this 
decay is almost absent when the kinetics of the total sig- 
nal is evaluated (Fig. 10(b)). This implies that the ESA 
contribution has a similar magnitude as the SE contribu- 
tion, and that they both decay as the population flows 
out of 64 and 65. This matches the observations of the 
experiments, where the kinetics of the low energy diag- 
onal peak show little dynamics even though it is clear 
from the pump-probe measurements that the dark states 
in the bottom of the manifold of excited states get pop- 
ulated on a 100 fs timescale. 

The changes in the shape of the cross-peak below the 
diagonal during relaxation from the upper to the lower 
band could provide valuable information on the relax- 
ation pathways. The experimental 2D spectrum at ^2 ~ 
45 fs (Fig. 2) shows that cp2 has its maximum at lower 
vi as compared to cpl. The simplest interpretation of 
this effect is that 65 gets populated from 67 while 64 gets 
populated from eg. The kinetics shown in Fig. 9 shows 
that 67 preferentially populates 64 and thus gives some 
support for such a conclusion. However, inspecting the 
simulated spectra in Fig. 2, we can conclude that the 
inclusion of ESA is essential to reproduce the tempo- 
ral evolution of the shape of the cross-peaks. The nega- 
tive ellipticity of the cross-peaks (an elongation along the 
anti-diagonal) remains after relaxation, but is only repro- 
duced when the ESA contribution is included in the simu- 
lations. The shape of the cross-peaks is thus not related 
to anti-correlated diagonal disorder'*® (site energies) or 
off-diagonal disorder (coupling disorder), but result from 
the overlap of the different signal contributions. 



IV. DISCUSSION 

A. Role of Charge Transfer States in [LuPc2] 

The ability to alter the CT state energies in [LuPc2] 
electrochemically gives a good opportunity to study the 
influence of CT states on the electronic structure and dy- 
namics in a coupled system. In this respect, it is instruc- 
tive to compare [LuPc2]^ to [LuPc2]* and for a moment 
neglect the influence of the counter-ion. It is already 
clear from Fig. 1 that the presence of the CT states in 
[LuPc2]^ dramatically changes the excited state struc- 
ture giving rise to new transitions and a redistribution 
of the oscillator strength. Comparing the evolution of 
the upper state(s) in [LuPc2]* and [LuPc2]^ illustrates 
that there is a strong impact on the dynamics as well: 
the lifetime decreases from 400 to 30 fs in the presence 
of the CT states. The results summarized in Fig. 8 show 
that our model can interpolate between [LuPc2]~ and 
[LuPc2]* by just shifting the energy of the CT states. 
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FIG. 9: Dynamics of eg (black solid), e^ (red dash), es (blue 
dot), 64 (magenta dash-dot), and ei (olive dash-dot-dot) after 
initial population of eg (a) and e^ (b). 
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FIG. 10: Kinetics of the simulated 2D spectra with (solid) 
and without (dashed) ESA contribution, (a) dpi (black) and 
dp2 (red), (b) cpl (black) and cp2 (red). 



and this allows us to pinpoint the mechanism for the 
slower dynamics in [LuPc2]*. Our model shows that the 
energetic splitting between the upper and lower band is 
significantly larger in [LuPc2]* (^ 2800 cm^^) as com- 
pared to [LuPca]" (-- 2000 cm^^). We find that the fast 
rates in [LuPc2]~ can only be reproduced when we in- 
clude vibrational (molecular) modes in the bath. From 
this we conclude that the 1600 cm^^ mode responsible for 
the vibrational sideband in [LuPc2] plays a central role 
in the relaxation processes. To fully explain the rates 
in the experiments we need to account for the influence 
of the counter-ion. The association with the counter- 
ion splits the degenerate energy levels and opens up ad- 
ditional relaxation pathways in the appropriate energy 
range [^ 1600 cm"-'^). Turning to [LuPc2]*; we find that 
the energy gap is too large for relaxation via the high fre- 
quency vibrational mode. Relaxation must thus proceed 
via the (weak) high energy wing of the Brownian oscil- 
lator spectral density used to describe the overdamped 
bath modes. The difference in the relaxation rates be- 
tween [LuPc2]^ and [LuPc2]* is thus not related to the 
enhanced coupling to the bath of the CT states, but de- 
termined by the differences in the energy level structure 
of the two systems. We note that the speedup of relax- 
ation in the presence of CT states depends on the details 
of the relaxation mechanism and the spectral densities. 
In the presence of high frequency vibrational modes, the 
larger system-bath coupling of the CT states may not 
necessarily lead to a significant increase in the relaxation 
rates as was shown recently in a simulation study of the 
photosystem II reaction center."*" 

From the geometrical structure and linear optical prop- 
erties, [LuPc2]~ seems to be a highly suitable candidate 
for investigations of fundamental processes like electronic 
coherence dynamics, dephasing and population relax- 
ation. Such studies would be valuable for a better under- 
standing of the more complicated dynamics observed in 
protein pigment complexes'*^"'*^ . However, the 2D spec- 
tra of [LuPc2]^TBA+ show distinct differences to the 
dimer models in the literature^^''*^''^'^"^^. These differ- 
ences cannot be explained by a rapid damping of the 
electronic coherences due to stronger system-bath cou- 
pling of the CT states. Rather, our analysis shows that 
including vibrational modes as well as employing a proper 
model for the two-exciton manifold is needed in order to 
understand the spectra and dynamics. The properties 
of the states in the two-exciton manifold responsible for 
ESA cannot be deduced from available quantum chemi- 
cal calculation with the needed precision, and therefore 
we need to resort to the exciton model. Perhaps the 
most prominent observation in the present experiments 
is the absence of a cross-peak below the diagonal at short 
population times (or lack of SE in pump-probe). While 
this requires a special structure of the two-exciton man- 
ifold, the same effect can show up in any coupled sys- 
tem. In fact, such cancelation happens readily in the 
exciton model when the two LEs have unequal transition 
dipole moments. The absence of the cross-peak should 
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thus not be seen as a direct effect of the CT states, but 
illustrates the importance of the two-exciton manifold 
and ESA for a proper interpretation of the spectra and 
dynamics. This might seem surprising because the 2D 
spectra are dominated by positive signal contributions 
(Fig. 2). However, ESA reshapes the peaks and strongly 
affects the observed kinetics (Fig. 10). We note that our 
simple model for the two-exciton band can account for 
the cancelation of the cross-peak below the diagonal as 
well as the absence of decay dynamics in the low energy 
band. To reach this agreement we need to shift the two- 
exciton band by —1000 cm~^ from the values estimated 
from the Frenkel exciton model. A similar redshift of the 
two-exciton manifold can be observed in [LuPc2]* (—550 
cm"-*^), indicating that electron correlation effects^^ are 
significant in this system even in the absence of the CT 
states. Inspecting the (full) pump-probe spectrum, one 
finds a very broad ESA covering the entire visible spec- 
tral range. ^^ In order to reproduce the ESA spectrum in 
the visible spectral range, we would need to include more 
(high energy) configurations to the two-exciton Hamilto- 
nian. The interaction of these high energy configurations 
with the ones included in our model could further con- 
tribute to the red-shift of the ESA observed in the exper- 
iments. 

The current experiments and analysis highlight the im- 
portance of the two-exciton manifold for the interpre- 
tation of time-resolved spectroscopy of processes which 
take place exclusively in the one-exciton manifold. In 
one-dimensional techniques, like pump-probe or transient 
grating, the presence of ESA may easily be overlooked 
leading to erroneous interpretation of the dynamics. In 
this respect, the combination of single- and double- 
quantum 2D spectroscopy is most useful to elucidate the 
role of ESA and the associated signal contributions. ^^'^* 

The pronounced influence of ESA in [LuPc2] makes 
in-depth analysis of the 2D spectra without numerical 
modeling difficult. Such simulations require a simple 
Hamiltonian mimicking the electronic structure of the 
system. There are many models which could fit the lin- 
ear spectrum and give rise to a suppressed cross-peak 
below the diagonal which grows as ^2 increases, e. g., 
the simple exciton model with unequal transition dipole 
moments of the two LEs or the isolated [LuPc2]^ model 
discussed above. However, none of these models are able 
to reproduce the dynamics in the cross-peak below the 
diagonal, or the double peak structure of the low energy 
band. By linking the exciton model to quantum chemical 
calculations, it is clear that precisely these features are 
signatures of special interactions, and that they provide 
important information on the dynamics in the system. 



B. Excited State Charge Transfer in 
[LuPca] TBA+ 

The association of [LuPc2]^ with the TBA+ counter- 
ion provides a natural explanation for the observation of 



the double peak structure in the low energy band in the 
2D spectra (Fig. 2 and Fig. 8). The formation of such ion- 
pairs has strong impact on the character of the excited 
states. The (positive) charge of the counter-ion serves to 
stabilize states representing CT towards the counter-ion. 
This gives rise to an asymmetry in the energies of the 
CT states, and as a result, the electronic eigenstates in 
[LuPc2]^TBA+ will refiect varying degrees of charge sep- 
aration. Relaxation between the different energy levels 
thus represents a net redistribution of electron and hole 
densities between the two rings, which can be followed 
by the evolution of the cross-peaks in the 2D spectra. 
By combining the analysis of the cross-peak dynamics 
within the exciton model with the results of the quan- 
tum chemical calculations, it becomes possible to obtain 
a molecular view of the CT dynamics. 

Fig. 11 shows the orbitals for the different bright tran- 
sitions with labels indicating the states and the corre- 
sponding features in the spectrum in Fig. 2. In addition, 
curved arrows are used to indicate the CT processes con- 
necting the different states. An excitation of the high 
energy band (dp3) populates an excited state with most 
of the charge density located on the side of the counter- 
ion (ring A). The system first relaxes via hole transfer 
(CT in the HOMO orbital) to ring B in 30 fs. This 
process is revealed via the faster rise of cp2 in the 2D 
spectrum. Further relaxation requires electron transfer 
(CT in the LUMO orbital) to ring B, and we find that 
this process also takes place on a 30 fs timescale. The 
system then goes to a stable (on ps timescale) charge 
separated state with a transition in the NIR. For this 
state, the excited state charge density is located on the 
ring closest to the counter-ion. Our pump-probe mea- 
surements reveal that population of the dark states in 
the NIR takes only about 100 fs. The combination of CT 
states and asymmetric interaction with the counter-ion 
thus opens up a relaxation pathway which very rapidly 
localizes the charge on one side of the complex. At first 
glance, one may assume that the counter-ion lowers the 
energy of both C^ and C^ and stabilizes the charge in 
both orbitals on the side of the counter-ion. However, 
it is clear from both the quantum chemical calculations 
as well as the exciton model that the molecular details 
of the interaction between [LuPc2]^ and the counter-ion 
need to be considered. The counter-ion induces different 
distortions of the molecular structure along the x- and y- 
coordinates, and this also affects the energies of the CT 
states representing CT into the corresponding orbitals on 
the macrocycle closest of the counter-ion. This selectivity 
is clearly manifested in the orbitals shown in Fig. 11, and 
shows that a detailed understanding of the ultrafast dy- 
namics requires careful consideration of the interactions 
in the ion-pair. 

[LuPc2]^TBA+ has many similarities to the so called 
special pair in the bacterial reaction center. The spe- 
cial pair is a dimer of bacteriochlorophylls coupled to 
CT states with higher transition energy than the funda- 
mental excitations. '^ However, because of different orien- 
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FIG. 11: Dominant orbitals involved in selected transitions for 
[LuPc2]"TBA+ using TD-BHLYP/SV(P). The legend to the 
left refers to the number of the state in Table I, and the legend 
to the right to the features assigned in Fig. 2. The timescales 
are obtained from the analysis of the kinetics of the 2D spec- 
tra. The curved arrows indicate the different electron- (ET) 
and hole-transfer (HT) processes taking place after excitation 
of the upper band in [LuPc2]~TBA+. 



tation of the transition dipoles, the lowest excited state 
(often named P^) carries most of the oscillator strength. 
Charge separation in the special pair takes place after di- 
rect excitation or energy transfer to P~ . Numerous works 
have concluded that P~ corresponds to a state with most 
of the charge density located on the M-side of the spe- 
cial pair. This localization of the charge density has 
been explained by specific interactions, such as hydrogen 
bonds^'*'^^, asymmetric interaction with the protein^ ^'^'^, 
or a special arrangement and interaction of the reaction 



center pigments 



20,59-61 



In the language of the exciton 



model, all these mechanisms serve to lower the CT states 
on one side of the complex in a similar way as the counter- 
ion in our case. Our results show that the presence of an 
asymmetric interaction facilitates relaxation to the bot- 
tom of the band, where the charge is localized on one 
side of the complex. In the special pair, the relaxation 
time between the upper and lower exciton band has been 



measured to be 65 fs.®^ The localization of the charge 
density on the M side of the special pair is thought to 
be important for the directionality of electron transfer 
in the bacterial reaction center, since the M side has a 
better overlap with the accessory bacteriochlorophyll on 
the L branch^^. Our experiments show that this type 
of charge localization can be mimicked by our model sys- 
tem, and followed in real time with 2D spectroscopy. The 
combination of 2D spectroscopy and theoretical analy- 
sis including density matrix propagation and quantum 
chemical calculations thus offers new possibilities to dis- 
entangle the complicated CT dynamics observed in the 
photosystem II reaction center^'^'^^. 



CONCLUSIONS 



In this work we have investigated the electronic struc- 
ture and excited state dynamics in a bis-phthalocyanine 
dimer ([LuPc2]~) where the resonance- and CT-coupling 
are of the same magnitude. To selectively investigate the 
effects from the resonance coupling, we have exploited 
the possibility to "turn off" the CT coupling by oxida- 
tion of the complex ([LuPc2]*). Based on linear absorp- 
tion, 2D spectroscopy and pump-probe measurements, 
we have shown that the mixing of exciton and CT states 
in [LuPc2] ~ leads to significant changes of the electronic 
structure and a tenfold speed-up of the excited state 
dynamics as compared to the excitonic [LuPc2]*. The 
speed-up of the dynamics in [LuPc2]~ can be qualita- 
tively understood from the changes in the energy level 
spacing due to the presence of CT states in combination 
with quasi-resonances between the energy gaps and high 
frequency vibrational modes. However, a detailed anal- 
ysis of the 2D spectra has shown that the interaction of 
[LuPc2] ^ with its TB A+ counter-ion needs to be consid- 
ered to account for all peaks and their evolution. Using 
quantum chemical calculations, we have demonstrated 
that the specific details of the interaction in the ion-pair 
determines the electronic structure of [LuPc2]~TBA+. 
The interaction with the counter-ion stabilizes CT states 
corresponding to CT into specific orbitals on the macro- 
cycle closest to the counter-ion. The subsequent im- 
balance in energy of the CT states results in an elec- 
tronic structure where the different excited states rep- 
resent varying degrees of charge separation. Relaxation 
in the excited state of the complex thus proceeds via 
electron- and hole-transfer processes, which give rise to 
distinct cross-peak dynamics in the 2D spectra. The 
analysis of the cross-peak dynamics, together with quan- 
tum chemical calculations, demonstrate the ability of 2D 
electronic spectroscopy to provide a detailed molecular 
view on these transient CT processes. 
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Appendix A. Quantum Chemical Calculations 

The BHLYP''''' functional was applied in Density Func- 
tional Theory (DFT) and time-dependent (TD-)DFT66 
calculations of optimal electronic ground state and low- 
est singlet excited-state geometries. This functional com- 
bines Becke's half-and-half exchange functional with the 
LYP correlation functional proposed by Lee, Yang and 
Parr. Among the tested functionals (PBE''^, BSLYP^^, 
and BHLYP^^), the overall best results were obtained us- 
ing BHLYP, which can be explained by the larger amount 
of exact exchange®^. The optical transitions were cal- 
culated using (TD)-DFT on the basis of the optimized 
geometries. We employed the split valence basis sets 
def-SV(P) (for N, C, H and O atoms) and def-ecp for 
lutetium'^° from the TURBOMOLE library. Herein, va- 
lence orbitals were calculated in a double-zeta basis and 
augmented with higher angular momentum polarization 
functions to account for the nonplanarity of the macrocy- 
cles. All quantum chemical calculations were performed 
with the Turbomole 5.7 package^^. 



Appendix B. Simulations of Excited State Dy- 
namics and Non-linear Spectra 

Excited state dynamics were simulated by apply- 
ing time-dependent Redfield theory to the diagonalized 
Frenkel exciton Hamiltonian of the system. ^^ Altogether 
eight singly excited electronic states were considered in 
the [LuPc2] dimer. In addition to the local Qx/Qy- 
excitations on each Pc ring, CT states corresponding to 
electron transfer from the excited state orbital on one 
macrocycle to the excited state orbital on the other ring 
were considered. The CT states were modeled as ex- 
cited states with zero transition dipole moment from the 
electronic ground state and increased reorganization en- 
ergy (with respect to the LE states). ^'^" The sign of the 
resonance coupling matrix elements, J ex, was determined 
from the structure. For the coupling between CT and LE 



states. Jet, we chose the signs of the matrix elements so 
that the Hamiltonian was invariant under rotation of the 
rings. The magnitudes of J^x and Jd were fitted to the 
linear absorption spectrum. Linear and non-linear spec- 
tra were calculated by a standard semiclassical response 
function theory. ^^ For calculation of non-linear spectra 
we constructed response functions corresponding to GSB, 
SE, and ESA in the Markov approximation^"^, i.e. ne- 
glecting correlations between the systems time evolution 
in different intervals of the response functions. Orienta- 
tional average was directly taken into account for each 
response function."*^ To describe line shapes and kinet- 
ics of the peaks, the dissipative population and coher- 
ence dynamics have to be calculated for the electronic 
states of the molecule. Dissipation was included via an 
analytical form of the time-dependent Redfield equation, 
invoking the secular approximation.^^ For the calcula- 
tion of the nonlinear spectra, the two-exciton manifold 
responsible for the ESA contributions needs to be in- 
cluded. As a first approximation we included the 28 
states which can be constructed by simultaneous exci- 
tation of two single exciton transitions (doubly excited 
states). This is a straightforward generalization of the 
usual Frenkel exciton model to the case of CT states. 
However, electron correlation effects^^, or inclusion of 
doubly excited monomeric states^'* can shift the energy 
of the two-exciton manifold. To account for these effects, 
not included in our model, we adjusted the energy of the 
two-exciton band to obtain the best agreement between 
simulation and experiment. A more rigorous treatment 
of the two-exciton CT band would have to account for 
various special properties of the combined doubly excited 
states, e.g., the fact that some double CT configurations 
result in normal double exciton states. ^'^ For calculation 
of the dephasing in the ESA part of the signal, the depop- 
ulation rates of the doubly-excited states are required. 
Here we assumed that the depopulation rates are equal 
in the one- and two-exciton manifolds. 

The coupling of the electronic transitions to the bath, 
described by the spectral density, determines both the 
line shapes and the population relaxation rates. Each 
transition was coupled to 3 vibrational modes with fre- 
quencies of 160, 725, and 1600 cm^^ and Huang-Rhys 
factors of S* = 0.3, 0.4, and 0.3, respectively. Addition- 
ally, one over-damped Brownian oscillator with reorgani- 
zation energy of 80 cm~^ and a decay rate of 100 fs~^ 
was used. For the CT states, we increased the coupling 
to the over-damped mode by a factor of 1.3. The vibra- 
tional modes are needed in order to reproduce the vibra- 
tional sideband in linear absorption as well as the fast 
rates found for [LuPc2]~. In all calculations of optical 
spectra, we explicitly averaged the signals over 1000 con- 
figurations where the transition energies were randomly 
sampled from a Gaussian distribution with a FWHM of 
300 cm-i. 
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I. SYNTHESIS 

A mixture of lutetium acetate (5 minol), 1,2-o-dicyanobenzene, and sodium carbonate 
in a stoichiometric ratio of 1 : 8 : 0.5 were stepwise heated to 280 °C and kept at con- 
stant temperature until sohdification of the meh signaled completion of the reaction. The 
reaction product was purified by sublimation under reduced pressure (< 1 Torr), following 
chromatographic separation and recrystallization.^'^ Column chromatography was carried 
out by dissolving the crude reaction mixture in a minimum amount of acetonitrile and fiush- 
ing it over aluminium oxide (AlO^c 90, neutral, activity grade I) using sodium methanolate 
(MeOH / MeO^Na+ (0.2 %)) as basic, polar mobile phase. The complex was precipitated 
by addition tetrabutylammonium bromide (TBABr) and crystallized for 3 days at -30°C. 
Recrystallization from DMF yielded [LuPc2]~TBA^ (5 %). All reagents were purchased 
from Sigma Aldrich and used as received without further purification. Solvents were applied 
in absolute forms. 

The complex was characterized using a number of spectroscopic techniques listed below. 

iR-NMR: (DMSO) 5 (ppm): 0.96 (t), 1.31 (m), 1.57 (m), 3.17 (s), 8.20 (q), 8.86 (q) 

i^C-NMR: (DMSO) 6 (ppm): 13.84 (CH3), 19.57 (CH2), 23.41 (CH2), 57.89 (CH2), 

121.91 (CH), 129.00 (CH), 137.15 (Q), 158.96 (Q) 

UV/VIS: (BN) Uraax (cm-i): 14245, 16180, 17700, 24330, 30030 

IR: (BN) Urnax (cm'i): 1112, 1330, 1381, 2875, 2935, 2965 

Ra: (BN) Umax (cm^i): 1505, 3050 

ESI-MS: (CH3OH/CH3CN) m/z: 1199.8 (M+) 



The radical was generated by dissolving [LuPc2]~TBA"'" in methanol and adding an equal 
amount of water and 0.5 %v/v concentrated hydrochloric acid. After evaporation of the 
solvent, the precipitate containing the neutral radical was isolated and dissolved in toluene. 

UV/VIS: (Toluene) u^ax (cm'^): 15190, 15950, 16830, 17630 
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FIG. 1: Spectrally resolved pump- probe signal (red circles) and the projection of the real part of 
the 2D spectrum (solid) for t2 = fs, 45 fs, and 400 fs (from bottom to top). The spectrum for 
t2 = 45 fs (400 fs) has been displaced by 2 (4) mOD along the vertical axis for clarity. 

II. PHASING OF THE 2D EXPERIMENT 

To extract the real part of the signal we employ a phasing procedure described previously.^ 
Accurate phasing is crucial for obtaining reliable information from the shape and dynamics of 
the various peaks. We find that all spectra need to be phased to their corresponding pump- 
probe spectrum to avoid errors due to the drift of the phase between the measurements. 
Fig. 1 shows the comparison of the pump-probe signal and the projection of the real part 
of the 2D spectrum for a few selected ^2 times. The phasing produces good agreement over 
the entire spectral range and for all times. 



HI. QUANTUM CHEMISTRY 

Normal mode analysis of [LuPc2]^ revealed a number of normal modes with significant 
Huang-Rhys factors with respect to the electronic state corresponding to the low energy 
band in linear absorption. The most striking mode is the 140 cm~^ (160 cm"^ in experi- 
ments) mode involving a modulation of the ring-to-ring distance. The motions of the most 
prominent normal modes are shown in Fig. 2 and the Huang-Rhys factors and frequencies 
are summarized in Table I. 



V cm ^ 


Huang-Rhys factor 


140 


0.5604 


726 


0.0697 


1573 


0.31 



TABLE I: Dominant ground-state vibrational frequencies and Huang-Rhys (HR) factors for 
[LuPc2]~ calculated with TD-DFT (BHLYP/SV(P)). 




140 cm-'' 



726 cm-1 



1573 cm-1 



FIG. 2: Motions of the atoms for the normal modes in Table I. 



Neither the exciton model nor the quantum chemical calculations are sufficient for a 
complete interpretation of the experimental results. The exciton model is needed for the 
interpretation of the excited state dynamics and simulations of the 2D spectra, but provides 
little insight into the molecular nature of the transitions. The quantum chemical calculation 
on the other hand provides a more realistic picture of the electronic structure of the molecule, 
especially in the presence of the counter-ion. Fig. 3 shows a comparison of the vertical 
transition energy for the first few excitations in the exciton model and the quantum chemical 
calculations with and without the counter-ion. The quantum chemical calculations have 
been shifted by —1200 cm^^. In the presence of a solvent, the two Pc rings are expected to 
planarize due to polarization effects leading to a red-shift of the transitions. The addition 
of the counter-ion leads to a splitting of the doubly degenerate bands and a red-shift of all 
transitions. We find that the splitting between the upper and lower band is somewhat larger 
in the calculations as compared to the exciton model. The similarities between the spectra 
of the two models suggest that 64 and 65 of the exciton model can be interpreted as ^3 and ^4 
from the quantum chemical calculations (as well as 67, es and gs, gg). A direct comparison of 
the transitions of the two models is complicated by the fact that the exciton model assumes 
that the ground state wave-function is evenly distributed over the two monomers. From the 
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FIG. 3: Transitions energies from the exciton model (blue) as compared to those of quantum 
chemistry for the [LuPc2]^TBA+ (red) and [LuPc2]~ (green). Both quantum chemical calculations 
have been shifted by —1200 cm^^. The labels show the respective states in the exciton model and 
quantum chemical calculations. Also shown is the linear absorption spectrum of [LuPc2]~. 

quantum chemical calculations, we can see that the interaction with the counter-ion leads to 
a significant localization of the ground charge density on the side of the counter-ion. Such 
effects are beyond our simplified model, and therefore we use the orbitals from the quantum 
chemical calculations to interpret the re- distribution of charge following relaxation. 

IV. TIME DEPENDENT REDFIELD RATES 

Relaxation processes in the first excited band of the dinier are described by the second 
order quantum master equation in the secular approximation. The system-bath interaction 
Hamiltonian is written in product form 



Hs-B = J2^'^nKn 



(1) 



where $„ (the energy gap operator) is the difference of the environmental potential energy 
surfaces of the ground state and the excited state \n). The operator if„ is the projector 
formed of the states \n), i.e. Kn = \n){n\. Starting with the well-known Nakajima-Zwanzig 
identity, the relaxation term in the equations of motion for the reduced density matrix p{t) 
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takes the form of an integral^ 

^ = -E r'"dr{Cr,niT)[K^,UsiT)K,At-r)U^ir)]_ 

Ot n "'0 

- Cnn{-r)[K^,Us{r)p{t - r)KnU^{r)]^} . (2) 

In the integral on the right-hand side we introduce slow envelop approximation Ug{t—T)p{t— 
T)Us{t — t) K. Ug{t)p(t)Us(t), by which the integro-difFerential Eq. 2 is converted into an 
ordinary differential equation with a time-dependent tensor. In the electronic eigenstate 
representation it reads 

^7 \relax = / , Rgbcdjt) Pcdjt) , (3) 

cd 

where 

Rahcd{t) = Kb E ^beedi^de] i) + hd 2_^ ^aeeci^ce] t) — ^cabdi'^db] t) — '^ dbaci'^ ca'-, t) ■ (4) 

e e 

The individual components of the tensor read 

Tabcd{uj-t) = Y.{a\Kn\h){c\Kn\d)Re /*"*" rfre^'^^C„„(r) . (5) 

„ Jo 



In addition, we neglect all terms coupling the evolution of coherences and populations (sec- 
ular approximation). In a long time limit, i.e. if to — ^ oo, this procedure leads to the 
well-known Redfield equations^. For our present purposes, we have to take into account the 
transient period, when the relaxation rates "build up", i.e. when they are time dependent. 
As we can see in Eq. 6, the Redfield tensor is composed of terms exhibiting the following 
form 

T{uj-t) = aRe f dTe-''^^C{T) , (6) 



'0 

where a is a constant and a; is a transition frequency in the excitonic band (we set to = 
0). To facilitate the calculations, we evaluate Eq. 6 for the case of a single overdamped 
Brownian oscillator with C{t) = C^e^"'^, where 7 is the inverse of the bath correlation time. 
We obtain 

,-^(l-e-'(cosM) + ^ 

7^ -|- cj^ ^ 7 

which leads in the long time limit to 



T{u- 1) = a^^ (1 - e-^*(cos(a;t) + - sin(a;t))) , (7) 



Too = a^ — ^ • 8 

7"^ -|- uj^ 



the Too are the components of the standard Redfield relaxation tensor. To capture the 
essential feature of the time dependence of the rates, we neglect the frequency of the tran- 
sition with respect to the bath correlation time, i.e. we assume cosut ^ l,sina;t ^ if 
g-7* > Q^ Neglecting the dependence of r{uj;t) on u, the Redfield tensor component is 
given by T{uj; t) = roo(l — e^'''*). For times longer than the bath correlation time we recover 
the standard (time independent) Redfield component Too- Assuming that the correlation 
functions of all components decay with the same constant 7, the time-dependent relaxation 
tensor reads 

R{t) = R^l - e-^') , (9) 

where Roo is the standard Redfield tensor. The form of the relaxation tensor implies that 
time dependent rates can combine fast relaxation (long time rates) with relatively narrow 
homogenous linewidths (initial rates). Furthermore, the ansatz in Eq. 9 greatly simplifies 
the calculation of the evolution superoperators needed to propagate the equation of motion. 
The evolution superoperator for the one-exciton band with the time dependent Redfield 
relaxation tensor has a form of a time ordered exponential 

U(t) = exp I - ^ / drRir)} . (10) 

*■ «, JO -" 

However, with our ansatz, Eq. 9 for the Redfield tensor, we can evaluate the propagator as 
an ordinary exponential 

U{t) = exp { - '-r J~'' + ^^-^ } , (11) 

which greatly speeds up the calculation. For the calculations of the rates we use a single 
overdamped Brownian oscillator with 7 = 0.01 fs~^. 

V. SUMMARY OF THE PUMP-PROBE RESULTS 

The pump-probe results for [LuPc2]"TBA''" presented in this study reveal similar results 
to those previously published by Prall et al. for the 13000-20000 cm^^ range^. Extending 
the probe range to cover the UV and NIR regions, ^'^ we can conclude that the dynamics in 
the B-band region is very similar to the Q-region. Just like in the previous study, we find 
that the spectra for t2 > 100 fs are independent of which band we pump. The transient 
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FIG. 4: Pump-probe spectra for [LuPc2]"TBA+ at ^2 = 100 fs (black solid), t2 = 500 fs (red 
dash), ^2 = 5 ps (blue dot), and ^2 = 15 ps (olive dash-dot). 
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FIG. 5: Decay associated difference spectra for [LuPc2] TBA+. The three time constants were 
3.8 ps (black solid), 15.8 ps (red dash), and 1 ns (blue dot). 

spectra in Fig. 4 show how the Q- and B-band regions (14000-16000 cm~^ and 30000 cm~^ 
respectively) GSB signals decay on a very similar timescale. 

We fitted the transient spectra to three exponentials and the decay associated difference 
spectra are shown in Fig. 5. The first component, with a 3.8 ps lifetime, has contributions 
from the GSB regions as well as ESA both in the NIR and VIS region. This also gives us 
the lifetime of the lowest state of the single exciton manifold. This excited state decays into 
a hot ground state, as evident by the dispersive lineshape of the second 15.8 ps component, 



representing cooling in the ground state. Similar conclusions were drawn by Prall et al.^, 
but can be re-enforced by the observation of the same type of features in both the Q- and 
-B band regions. An additional very slow component was needed to account for residual 
signals. However, as can be seen from Fig. 5, the amplitude of this component is very low. 
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